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. Summary. In the first part of this chapter we analyze the contact intensity C, 

{3jQf which has been introduced by Tan [Ann. Phys. 323, 2952 (2008)] and appears in 

several physical observables of the strongly correlated two-component Fermi gas. 

■ We calculate the contact C in the full BCS-BEC crossover for a uniform superfluid 

■ Fermi gas by using an efflcient parametetrization of the ground-state energy. In the 
^ , case of harmonic confinement, within the Thomas- Fermi approximation, we derive 

analytical formulas of C in the three relevant limits of the crossover. In the second 
part of this chapter we discuss the extended superfluid hydrodynamics we have 
recently proposed to describe static and dynamical collective properties of the Fermi 
gas in the BCS-BEC crossover. In particular we show the relation with the effective 
theory for the Goldstone fleld derived by Son and Wingate [Ann. Phys. 321, 197 
(2006)] on the basis of conformal invariance. By using our equations of extended 
^ ■ hydrodynamics we determine nonlinear sound waves, static response function and 

, structure factor of a generic superfluid at zero temperature. 



>■ ■ 1 Contact intensity 

■ It has been proved by Tan that the momentum distribution pa {k) in an arbitrary 

system consisting of fermions in two spin states {a =t, i) with a large scattering 
length has a tail that falls off as 



for fc CO, where C is the so-called contact intensity IJ. Here large s-wave scattering 
length a means that |a| 3> tq, where ro is the effective interaction radius. Under this 
condition Tan '1[ has shown that C is related to the total energy E of the Fermi 
system by the rigorous expression 

^^W^ (2) 
da ' ^ ^ 



where the derivative is taken under constant entropy and, in general, C depends 
on the number A'' of fermions, the scattering length a and the parameters of the 
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trapping potential [21 El ■ Remarkably, Eqs. ((T]) and ((21 work also at finite temper- 
ature and in this case C will be a function of T [2l |4]. Tan has also derived, for 
finite scattering lengths, a generalized virial theorem and a generalized pressure re- 
lation where the contact C appears [3]. The contact intensity C appears also in 
other physical observables of the strongly correlated Fermi system. For instance, the 
radio- frequency spectroscopy shift is proportional to C ^Sj [6] [T] [8], and the same 
happens to the photoassociation rate [5]. Very recently, it has been shown that the 
contact C gives the asymptotic tail behavior of the shear viscosity as a function of 
the frequency [lO] . 

Using the methods of quantum field theory, Braaten and Platter have rederived 
[111 112] the Tan's universal relations [U (2] [3] . In addition, they have shown that the 
contact intensity can be written as 

j /(V't+(r)V;;+(r)4(r)^Hr)) d\ , (3) 

where ■(/'^(r) is the fermionic field operator of spin a and g = Ana/il — {2akcut / 1^)) 
is the bare coupling constant of the Fermi pseudo-potential interaction, with kcut 
the ultraviolet wavenumber cutoff [111 I13| . Braaten and Platter have also shown 
that the number Npairijc) of pairs of fermions with opposite spins in a small ball 
of volume 47rs^/3 centered at the point r scales as s''C(r)/2 for s — >■ 0, where 
C(r) = g^{7/)+(r)7/)+(r)V;;(r)i/)-t.(r)) is the contact density [TT1IT3] . 

Explicit expressions of the universal quantity C have been derived by Tan [1] 
for a uniform superfluid Fermi gas at zero temperature only in three limits: the 
Bardeen-Cooper-Schrieffer (BCS) limit of weakly bound Cooper pairs, the unitarity 
limit of infinite scattering length, and the Bose-Einstein condensate (BEC) limit of 
weakly-interacting bosonic molecular pairs. 

In this section we calculate the contact C as a function of the inverse scattering 
parameter l/ikpa) for a uniform superfiuid Fermi gas in the full BCS-BEC crossover. 
To perform this calculation we use an efficient fitting formula of the ground-state 
energy [251 126[ and the Tan's equation ((2]). We find that the contact C has a max;- 
imum close to the unitarity limit of infinite scattering length. We also consider the 
interacting Fermi system under harmonic confinement. For this superfluid Fermi 
cloud we derive analytical formulas of the contact C in the three relevant limits of 
the crossover. 

1.1 Uniform superfluid Fermi gas at zero temperature 

In the case of a zero-temperature uniform two-component superfiuid Fermi gas of 
total density n = n-^ + {n^ — 714,), large scattering length a (a 3> ro) in a volume 
V, the energy density can be written as 

V = h'pfii^) . (4) 

V 5 Kpa 

where f{y) is a universal function of inverse interaction parameter y — l/ikpo), with 
tF = h^kp/{2m) the Fermi energy and kp ~ (STv^n)^^'* the Fermi wave number[4]. 
We observe that at finite finite temperature T the function f(y) is substituted by a 
more general universal function (l>{y,t), where t = T/ep is the scaled temperature. 
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$(y, t) has been studied with Monte Carlo methods by Bulgac, Drut, and Migierski 
[141 115) . but only in the unitarity limit [y = Q). 

It is straightforward to derive from Eq. ([2]) the expression of the contact density 

The behavior of f{y) is well known in three relevant limits: 

ri + f i + 0(l/y') , 
f{y)={ ?-Cy + o(2/'), (6) 

I W-.\+0{l/y''^), y»l 

In fact, in the weakly attractive limit {y ^ —1) one expects a BCS Fermi gas of 
weakly bound Cooper pairs. As the superfluid gap correction is exponentially small, 
the function f{y) follows the Fermi-gas expansion |16II17| . In the so-called unitarity 
limit {y = 0) one expects that the energy per particle is proportional to that of a 
non-interacting Fermi gas with the universal constant ^ given by ^ ~ 0.42 [18) . Note 
that more recent auxiliary-field Monte Carlo results fT5] predict a smaller value for ^, 
namely ^ ~ 0.38, while the experiment performed at Ecole Normale Superieure [20| 
suggests ^ ~ 0.40. The first correction to this behavior, shown in Eq. (O, has been 
estimated from Monte Carlo data with ~ 1 [21[ . In the weak BEC limit {y ^ 1) 
one expects a weakly repulsive Bose gas of dimers. Such Bose-condensed molecules 
of mass niM = 2m and density um = n/2 interact with a positive scattering length 
fflM = Va with V — 0.6, as predicted by Petrov et al. [25]. In this regime, after 
subtraction of the molecular binding energy, the function f(y) follows the Bose-gas 
expansion [23| . It is easy to obtain the contact intensity C by using Eqs. © and ([5} 
in the relevant limits of the crossover. One finds 

C |fc|,na^ -f O(a^) , y<-l 
= ) ^kpnC + 0{l/a) , -1 « y « 1 (7) 
^ [|4na2+0(a^/^), y » 1 

in agreement with the previous determinations of Tan yj and Werner, Tarruell and 
Castin [9]. Notice that in the BEC limit we have removed the binding energy of 
molecules. Moreover, very recently finite temperature corrections to Eq. (O were 
given in Ref. [24]. 



1.2 Contact intensity in the BCS-BEC crossover 

Now we want to calculate the behavior of C in the full BCS-BEC crossover. In 2005 
we have proposed [25] the following analytical fitting formula 

f{y) = ai - 02 arctan (03 y iT-TTl] ' 

interpolating the Monte Carlo energy per particle [18| and the limiting behaviors 
for large and small \y\. Eq. ^ is very reliable [251 126] and it has been successfully 
used by various authors for studying grund-state and collective properties of this 
superfluid Fermi system [271 1281 [29] 1301131] . The parameter ai is fixed by the value 
^ of f(y) at y = 0, the parameter 02 is fixed by the value of f{y) at y — 00, and 
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03 is fixed by the asymptotic 1/y coefficient of e{y) at large \y\. The ratio /3i//32 
is determined by the linear behavior C, of e(y) near 3/ = 0. The value of /3i is then 
determined by minimizing the mean square deviation from the Monte-Carlo data. 
Of course, we have considered two different set of parameters: one set in the BCS 
region {y < 0) and a separate set in the BEC region {y > 0) ^251. Table 1 of [25] 
reports the values of these parameters, with = = 1 in the BCS region but 
C = C+ = 1/3 in the BEC region. 
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Fig. 1. Upper panel: Scaled ground-state energy E/CVep) of the uniform Fermi 
gas as a function of the inverse interaction parameter l/{kFa) in the BCS-BEC 
crossover. Lower panel: Scaled contact intensity C / (Vkpn) of the Fermi system as 
a function of the inverse interaction parameter l/{kFa)- Here V and n = N/V are 
the volume and the density of the Fermi gas, kp — (STv'^n)^^^ is the Fermi wave 
number and ep ~ h^kp/{2m) is the Fermi energy. Two different parametrization of 
the universal function f{y): solid lines are obtained with Eq. ((8|, while dashed lines 
are calculated using the Fade approximant of [52] • 

Here we use Eq. ([8]) to calculate the contact density given by Eq. ([5|, but contrary 
to [25], we choose = = 1 to ensure the continuity of f'{y) at y = 0. 

Notice that the recent experimental results obtained at Ecole Normale Superieure 
[20j indeed suggest the continuity of f'{y) at j/ = 0. In this way, in the BEC region 
(32 = 0.1517 while /3i is unchanged (see Table. 1 of 25j). In the upper panel of Fig. 
1 we plot the ground-state energy E while in the lower panel we plot the contact C, 
both as a function of the inverse interaction parameter y — l/ikpa). For comparison, 
in addition to the data obtained with our method (solid lines), we insert also the 
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results (dashed lines) one obtains using the Pade parametrization of f{y) proposed 
by Kim and Zubarev [32]. The figure shows that solid and dashed lines are always 
close each other, apart for — 1 <C y < where our fitting formula is smoother (and 
closer to the Monte Carlo data 18 ). Moreover, the scaled contact C / (Vkpn) as 
a function of l/{kFa) has its maximum near to the unitarity limit l/(kFa) = 0: 
the position of the maximum is located at l/(kFa) ~ 0.05. Remarkably, the scaled 
contact has a behavior quite similar to the Landau's critical velocity Vc (at which 
there is the breaking of superfluid motion), calculated along the BCS-BEC crossover. 
In fact, also «c goes to zero for y — > ±oo and it has a peak at j/ ~ 0.08 |33j . Clearly, 
the contact C exhibits a maximum close to unitarity because we have subtracted 
the molecular binding energy contribution, given by in the BEC side (j/ > 0). 
Including this energy term it is easy to show that C increases monotonically from 
the BCS side to the BEC side, according the definition in Eq. Nevertheless, as 
prevously stated, the radio-frequency spectroscopy shift is proportional to C, and 
its maximum around unitarity has been shown in Ref. 6 , where the same artificial 
subtraction has been adopted. 

1.3 Trapped superfluid Fermi gas 

Let us now consider the superfluid Fermi gas under an external harmonic confine- 
ment 

f/(r) = ima;"(a;' + y' + ^') . (9) 

In the limit of a large number N of fermions we can use the local density (Thomas- 
Fermi) approximation [Ml [33 ESI 133 [Ml [33 HOI SI] and the energy of the system 
can be written as 

E = j ||n(r)e^(r)/(-^-l^)+n(r)C/(r)} d\ , (10) 

where ^f{v) = ft^fcF(r)^/(2m) is the local Fermi energy and kF = (37r^n(r))^''^ is 
the local Fermi wave number. As in the uniform case, this expression is not very 
useful without the knowledge of the universal function f(y). 

The numerical calculation of the contact indensity C for a harmonically trapped 
Fermi superfluid in the full BCS-BEC crossover by using Eqs. ([2]) and ([8)l is very 
demanding. Consequently, we calculate C only in the three relevant limits of the 
BEC-BEC crossover. In these limits we obtain useful analytical expressions for the 
contact C. 

BCS limit. In the BCS hmit (a 0") from Eqs. © and ^ we find 

f = fS/^-W"(r)^-W>^'^' 

where n(r) is the density profile of the ideal Fermi gas in the harmonic potential 
©, given by [42] 

where tf = V2{iNY'^' an is the Fermi radius of the cloud, with an ~ \/'h/{mLu) 
the characteristic harmonic length. Inserting this density profile into Eq. pif) and 
using Eq. ([SJ we get the contact intensity in the BCS limit: 
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o-'^^(-ymf'- (13) 

283571 an \aH J 

Unitarity limit. In the unitarity limit (a — >■ ±oo) from Eqs. ((Gjl and (|10|l we have 
where n(r) is the density profile of the unitary Fermi gas in the potential ([9]), namely 

EH 



where tf = ^/25^''^(3iV)^'"' aj/ is the Fermi radius of the unitary cloud. Inserting 
this density profile into Eq. (|14p and using Eq. ^ we obtain the contact in the 
unitarity limit: 

BEC limit. In the BEC limit (a — >■ 0"'") it is straightforward to calculate the contact 
C. In fact, the explicit formula of the ground-state energy of the dilute BEC is well 
known [43] and for a BEC of molecules it is given by 

where V ~ 0.6 and N/2 is the number of molecules. Then from Eq. ^ the contact 
intensity reads 

27r^ /iSpy/"^ (oNV''' 



„2 



2 Extended superfluid hydrodynamics 

In this section we discuss the extended Lagrangian density of superfiuids we have 
proposed few years ago [31] and applied to study mainly the unitary Fermi gas [341 
13511361 [37l 1381 l39l 1401 |4T] . The internal energy density of this Lagrangian contains a 
term proportional to the kinetic energy of a uniform non interacting gas of fermions, 
plus a gradient correction of the von-Weizsacker form A?i^/(8m)(Vn/n)^ j44|. This 
approach has been adopted for studying the quantum hydrodynamics of electrons by 
March and Tosi [15], and by Zaremba and Tso [IB]. In the context of the BCS-BEC 
crossover, the gradient term is quite standard [25l l32]|47ll48ll49l[50l[5T][52l[53] . 
In particular we show the relation between our approach and the effective theory 
for the Goldstone field derived by Son and Wingate [51], and improved by Manes 
and Valle [55], on the basis of conformal invariance. Finally, by using our equations 
of extended superfluid hydrodymics at zero temperature we calculate sound waves, 
static response function and structure factor of a generic superfluid. 
The extended Lagrangian density of superfiuids is given by 



t2 

C = -nen (Vef n-£{n, Vn) - U{r) n , 



(19) 
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where n{r,t) is the local density, U{r) is the external potential acting on particles, 
and m the mass of superfluid particles. In the case of superfluid bosons 0{r, t) is the 
phase of the condensate order parameter, while in the case of superfluid fermions 
0{Y,t) is half of the phase of the condensate order parameter (of Cooper pairs). 
£{n,'Vn) is the internal energy density of the system. Note that we are supposing 
that this equation of state £{n, Vn) can depend not only on the local density n(r, t) 
but also on its space derivatives. For this reason we call p9p the extended superfluid 
Lagrangian. We stress that in the context of the BCS-BEC crossover the extended 
internal energy density could be written as 

£{n,Vn)=£o{n)+\^^^^^ , (20) 
om n 

where 

£oin) = |n.^/(-i-) (21) 

is the energy density discussed in the previous section (see Eq. (Q) which depends 
on the universal function f(y) of the BCS-BEC crossover, while the second term is 
the gradient correction of the von-Weizsacker form [35]. In the BCS-BEC crossover 
we expect that 1/6 < A < 1/4, where A = 1/6 is the appropriate value in BCS 
regime of weakly-interacting superfluid fermions of mass m [451 146| , while A = 1/4 
is the appropriate value in the deep BEC regime of weakly-interacting superfluid 
bosonic dimers of mass 2m [35| . 

By using the Lagrangian density (|19[) the Euler-Lagrange equation 

^_1^_V.^^=0 (22) 

gives 

^ + -V ■ (n Ve) = . (23) 
at m 

The Euler-Lagrange equation 

dC d dC dC 
g£.__ga^. PL. ^ 

dn dtdh 9(Vn) ^ ' 



gives instead 

2m 



he+^{Vef + U{r) + X{n,Vn)=Q , (25) 



where 

X(n,Vn)^f -V.^ (26) 

is the local chemical potential of the system (see also [53] and [55)). 

The local field velocity v(r, t) of the superfluid is related to the phase 9(r,t) of 
the condensate by 

v(r,t) = -Ve(r,t) . (27) 
m 

This definition ensures that the velocity is irrotational, i.e. V A v = 0. By using the 
definition (|27[) in both Eqs. (|23|) and (|25|l and applying the gradient operator V to 
Eq. ([25[) one finds the extended hydrodynamic equations of superfluids 
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^ + V ■ (n v) = . (28) 



9v „ 



1 ■> 

-mv + U{r) + X{n, Vn) 



. (29) 



We stress that in the presence of an external confinement U (r) the chemical po- 
tential jj, of the system does not coincide with the local chemical potential X{n, Vn). 
The chemical potential /i can be obtained from Eq. (I25p setting 9(r,t) — —fit/h and 
v(r, t) = 0, such that 

U{r) + X{no,Vno)= , (30) 

where no(r) is the ground-state local density. 

The Lagrangian density (|19p depends on the dynamical variables 6{r, t) and 
n(r, t). The conjugate momenta of these dynamical variables are then given by 

d£ 

TTg — r — ~hn , (31) 

09 

(32) 

and the corresponding Hamiltonian density reads 

T-L — -Ke 9 + Tinfi — C = —h n9 — C , (33) 

namely 

■H^—{V9fn + £{n,Vn) + U{v)n, (34) 
2m 

which is the sum of the flow kinetic energy density ?i^(VS)^n/(2m) = [l/2)mv^n, 
the internal energy density £{n, Vn), and the external energy density U{v)n. 



2.1 Extended hydrodynamics in terms of Goldstone field 

Note that taking into account Eq. (|25p one immediately finds 

X{n, Vn) n^-h9n- —{V9f n - U{v) n . (35) 
2m 

Consequently the Lagrangian density (|19|) can be rewritten as 

C = X{n,Vn)n~ £{n,Vn) . (36) 

Remarkably 

P{n,Vn) = X{n,Vn)n- £{n,Vn) (37) 

is the local pressure of the system as a function of the density and its spatial deriva- 
tives, which can be written as a function of the local chemical potential X and its 
spatial derivatives, namely 

C = P(X, VX) . (38) 

This result, based on a Legendre transformation, is clearly illustrated in the book of 
Popov [SS] and used in the recent papers of Son and Wingate 54 and Manes and 
Valle [55] ■ Finally one can introduce the Goldstone field (/>(r, t) as 

^{v,t)^9{v,t) + ^t . (39) 
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In this way, by using again Eq. (I25p . one can write 

X^fi-h^-^{V^f-U{r), (40) 

Thus, the Lagrangian density (|38[) actuaUy depends only on the Goldstone field 
(j}{r,t). This is exactly the main message of the paper of Son and Wingate [54], 
which however traces back to the older results of Popov 56 . 



2.2 Application: the Unitary Fermi gas 

Let us now suppose that the equation of state of the superfluid at zero temperature 
is that of the unitary Fermi gas, i.e. 

£ir^^^n) = - — {2,.)'^n' +X——-, (41) 
where ^ ==; 0.4 and A ~ 0.25 [371 E]. It follows from Eq. ^ that 

x(„,v„).|;P,")"'«„--a|1<X2)!, ,«) 

by taking into account that the surface terms give zero contribution. In addition we 
get from Eq. ((33 that 

r,^ TV \ 2 2N2/3f 5/3 , (Vw)^ , , 

P{n,'Vn) ^ - — {Sn ) ' ' - X—- -, (43) 

5 2m 4m n 



The Lagrangian density of Eq. (I38|l is then obtained by finding n and Vn as functions 
of X and VX by inverting Eq. (|42|l . This can be done in terms of a derivative 
expansion. One gets n = (2m^)^/^X^''^/(37r^R^) and 

C^Clo+ Cnlo (44) 

where 

Clo=co^X'^\ (45) 

with 

25/2 

^0 = ' (46) 

is the Lagrangian density at the leading order, and 

n^(yx_ 
h y/x 



ml/2 (y^)2 

jCnlo = Cl -— ^ ' , (47) 



with 

3 • 2^/2 

^1 - -^^^ ' (48) 

is the next-to-leading contribution to the Lagrangian density. The Lagrangian den- 
sity (|44p is the same of that derived by Son and Wingate [S3] from general coordinate 
invariance and conformal invariance. Actually, at the lext-to-leading order Son and 
Wingate have found an additional term [54], which has been questioned by Manes 
and Valle [5S] and is absent in our approach. 
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2.3 Nonlinear sound waves, static response function and structure 
factor 

In this subsection we consider the following zero-temperature equation of state of a 
generic superfiuid 

£{n,Vn)=£o{n) + \-^^^j^ . (49) 

Here the internal energy is the sum of two contributions: a generic internal energy 
Soin) which depends only of the local density n{r,t) (for instance that of Eq. (|21ll ') 
plus the gradient correction of the von Weizsacker type, where the coefficient A can 
be a function of the interaction strength. 
The equation of motion H29|) becomes 



5v „ 



imv" + U{r) + Xo{n) ~ X^ — 
A Am V 



= , (50) 



with 

Mn) ^ f . (51) 

We are interested on the propagation of sound waves in superfluids. For simplicity 
we set 

U{r) = , (52) 

and consider a small perturbation h{r, t) around a uniform and constant configura- 
tion no, namely 

n(r, t) = no + n(r, t) . (53) 

Neglecting quadratic terms in n and v we derive the linearized equations of extended 
superfiuid hydrodynamics 

^ + noV • V = , (54) 
ot 

no|^+c?Vn-A^V(V2n)=0, (55) 

where Cs is the sound velocity, given by 

^2 _ no dXojno) _ np d'^£o{no) 
" m dn m dn^ 

Applying the operator ^ to Eq. (|54p and the operator V to Eq. (|55|l and subtracting 
the two resulting equations we get 

&2-"sV'+\^,V')n{r,t) = Q. (57) 

This is the wave equation of the small perturbation h{Y, t) around the uniform 
density no. We stress that the effect of the gradient term in the equation of state is 
the presence of quartic spatial derivatives in this wave equation. 

It is straightforward to show that the wave equation admits the real solution 

n(r, t) = A e^t'i--"''*) + A* e-'(i—"5*) , (58) 
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where the frequency u)q and the wave vector q are related by the Bogoliubov-like 
dispersion formula 

or equivalently 



1^9 = Ca«7\/l + aq^ = Csq \^ + -q^ - — g"* + ...j , (60) 

with a — \h? /{Am?cl). Thus, the dispersion relation u)q is linear in g = |q| only for 
small values of the wavenumber q and becomes quadratic for large values of q. 

We observe that, for a generic many-body system, the dispersion relation can 
be written as [ST] 

where mn{q) is the n moment of the dynamic structure function S{q,u}) of the 
many-body system under investigation, i.e. 

mn{q) — du! S{q,Lu) (hcu)" . (62) 
Jo 

Note that Eq. H61|) is not exact and is valid under the approximation of a single- 
mode density excitation. It therefore only gives an upper bound for the dispersion 
relation. This is important for large q, for which quasiparticles other than density 
excitations may contribute. In our problem we have 

= ^ (63) 

and 

In general, the static response function x{q) is defined as [57] 

X(g) = -2m_i(g) , (65) 

and in our problem it reads 

2 

x(<?) = TT^ . (66) 

2m ^ 

or equivalently 

Xil) = -^TT-^ = (1 - + + ...) , (67) 

mci 1 + aq^ mci ^ ' 

where again a = \h^/{4m^c^). 

The static structure factor S{q) is instead defined as |57] 

poo 

Siq) = mo (a) = duj S{q, uj) , (68) 
Jo 
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but it can be approximated by the expression 



S{q) = y^mi{q)m-i{q) , 



(69) 



which gives an upper bound of S{q), i.e. 5(g) > S{q) 57 . In our problem we 
immediately find 



Our results clearly indicate that it should be possible to observe experimentally the 
effect of the dispersive von-Weizsacker-like gradient term from sound-wave measure- 
ments. 



3 Conclusions 

In the first part of this contribution we have calculated the contact C as a function 
of the inverse scattering parameter l/(fcFa) for a uniform superfluid Fermi gas in 
the full BCS-BEC crossover at zero temperature. We have found that the contact C 
has a maximum close to the unitarity limit of infinite scattering length, in analogy 
with the behavior of the Landau's critical velocity Vc, at which there is the breaking 
of superfluid motion [33] . We have also considered the interacting Fermi system un- 
der harmonic confinement. In this case, we have derived analytical formulas of the 
contact intensity C in the three relevant limits of the crossover. Our results can be 
experimentally tested with ultracold atomic clouds by measuring one of the quanti- 
ties which are directly related to the contact intensity C: the tail of the momentum 
distribution, the derivative of the total energy with respect to the scattering length, 
the radio-frequency spectroscopy shift, or the photoassociation rate. In the second 
part we have analyzed some properties of the extended superfluid hydrodynamics 
[34) : in particular, we have shown its strict relation with the low-energy effective fleld 
theory built on the Goldstone mode. Finally, by using the extended hydrodynamics 
we have calculated, for generic superfluid in the absence of external confinement, 
the nonlinear dispersion relation of sound waves, and, as a by-product, both static 
response function and structure factor. 
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